MS221 CB B 


pen 
rsi 


TheO 
Unive 


A second level 


interd 


e 
g 
2 


: 
: 


isc 


MS221 CBB 


2 


ee) 


__ Asecond level 
inferdisciplinary —s_ 
course 


BLOCK B 
EXPLORING ITERATION 


Computer Book B 


ICS 


oring 


Maeathemat 


Exp! 


Prepared by the course team 


About this course 


This computer book forms part of the course MS221 Exploring Mathematics. 
This course and the courses MU120 Open Mathematics and MST121 Using 
Mathematics provide a flexible means of entry to university-level mathematics. 
Further details may be obtained from the address below. 


MS221 uses the software package Mathcad (MathSoft, Inc.) to investigate 
mathematical concepts and as a tool in problem solving. This software is 
provided as part of the course. 


This publication forms part of an Open University course. Details of this and 
other Open University courses can be obtained from the Course Information and 
Advice Centre, PO Box 724, The Open University, Milton Keynes, MK7 625, 
United Kingdom: tel. +44 (0)1908 653231, e-mail general-enquiries@open.ac.uk 


Alternatively, you may visit the Open University website at 
http: //www.open.ac.uk where you can learn more about the wide range of 
courses and packs offered at all levels by The Open University. 


To purchase a selection of Open University course materials, visit the webshop 
at www.ouw.co.uk, or contact Open University Worldwide, Michael Young 
Building, Walton Hall, Milton Keynes, MK7 6AA, United Kingdom, for a 
brochure: tel. +44 (0)1908 858785, fax +44 (0)1908 858787, e-mail 


ouwenq@open.ac.uk 


The Open University, Walton Hall, Milton Keynes, MK7 6AA. 
First published 1997. Second edition 2002. New edition 2004. 
Copyright © 2004 The Open University 


All rights reserved; no part of this publication may be reproduced, stored in a 
retrieval system, transmitted or utilised in any form or by any means, electronic, 
mechanical, photocopying, recording or otherwise, without written permission from 
the publisher or a licence from the Copyright Licensing Agency Ltd. Details of such 
licences (for reprographic reproduction) may be obtained from the Copyright 
Licensing Agency Ltd, 90 Tottenham Court Road, London W1T 4LP. 


Open University course materials may also be made available in electronic formats 
for use by students of the University. All rights, including copyright and related 
rights and database rights, in electronic course materials and their contents are 
owned by or licensed to The Open University, or otherwise used by The Open 
University as permitted by applicable law. 


In using electronic course materials and their contents you agree that your use will 
be solely for the purposes of following an Open University course of study or 
otherwise as licensed by The Open University or its assigns. 


Except as permitted above you undertake not to copy, store in any medium 
(including electronic storage or use in a website), distribute, transmit or re-transmit, 
broadcast, modify or show in public such electronic materials in whole or in part 
without the prior written consent of The Open University or in accordance with the 
Copyright, Designs and Patents Act 1988. 


Edited, designed and typeset by The Open University, using the Open University 
TRX System. 


Printed and bound in the United Kingdom by The Charlesworth Group, 
Huddersfield. 


ISBN 0 7492 6648 1 
3.1 


Contents 


Guidance notes 


Chapter B1 
Section 4 Iterating real functions with the computer 
4.1 Sequences obtained by iterating f(x) = 2* +c 
4.2 The bifurcation diagram 
4.3 Other quadratic iteration sequences (Optional) 
4.4 Sequences obtained by iterating real functions 


Chapter B2 
Section 5 Visualising affine transformations 
5.1 Linear transformations 
5.2 Affine transformations 


Chapter B3 
Section 5 Iterating linear transformations with the 
computer 


5.1 Iterating linear transformations 
5.2 Iterating affine transformations 
5.3 Iterated function systems (Optional) 


Solutions to Activities 


or ot Oot 


11 
14 
15 


17 
17 
17 
te 


24 


24 
24 
30 
31 


34 


Guidance notes 


This computer book contains those sections of the chapters in Block B 
which require you to use Mathcad. Each of these chapters contains 
instructions as to when you should first refer to particular material in this 
computer book, so you are advised not to work on the activities here until 
you have reached the appropriate points in the chapters. 


In order to use this computer book, you will need the following Mathcad 


files. 


Chapter B1 

221B1-01 Iterations of f(x) =z?+c¢ 

221B1-02 Overview of iterations of f(x) = x? +c 
221B1-03 Iterating real functions 


Chapter B2 


221B2-01 Linear transformations 
221B2-02 Affine transformations 


Chapter B3 

221B3-01 Iterating linear transformations 
221B3-02 Iterating affine transformations 
221B3-03 Iterated function systems (Optional) 


Instructions for installing these files onto your computer’s hard disk, and 
for opening them, are given in Chapter AO of MST121. 


The computer activities for each chapter also require you to work with 
Mathcad documents which you have created yourself. 


Activities based on software vary both in nature and in length. Sometimes 
the instructions for an activity appear only in the computer book; in other 
cases, instructions are given in the computer book and on screen. 
Feedback on an activity is sometimes provided on screen and sometimes 
given in the computer book. 


For advice on how each computer session fits into suggested study 
patterns, refer to the Study guides in the chapters. 


Chapter B1, Section 4 
Iterating real functions with the computer 


In this section you will explore iteration sequences, using Mathcad 
worksheets that draw graphical iteration diagrams and calculate related 
information, such as fixed points and gradients. 


In Subsection 4.1, you will explore sequences obtained by iterating 

quadratic functions of the form f(x) = x* +c. You will see that many 

different types of behaviour occur; for example, such sequences can tend to 

p-cycles for various periods p. This subsection is quite long; if you are The concept of a p-cycle with 


short of time, then concentrate on Activities 4.1 and 4.2. p > 2 was defined in the main 


text. 
In Subsection 4.2, Mathcad is used to provide an overview of how the 


behaviour of these iteration sequences changes as c varies, and you will see 
that the changes in behaviour are surprisingly complicated. Subsection 4.3 
is optional; it contains a discussion of certain other quadratic iteration 
sequences. In Subsection 4.4, you will use a Mathcad worksheet that allows 
you to study the sequences obtained by iterating any real function f. 


4.1 Sequences obtained by iterating f(x) = x’ +c 


In this subsection you will explore sequences obtained by iterating 
functions of the form f(x) = x? +c, where c is a parameter whose value 
can be varied. 


You have seen that if we wish to study the long-term behaviour of 
iteration sequences, then it is useful to find and classify any fixed points 
and 2-cycles of the function. 


The fixed point equation and the 2-cycle equation of the function 

f(z) =x? +c can be solved to give formulas, in terms of c, for the fixed 
points and for those points that are members of a 2-cycle. The Mathcad 
worksheet that accompanies this subsection uses these formulas to 
calculate any fixed points and 2-cycles of f. The number of fixed points 
and 2-cycles depends on the value of c, but f always has at most two fixed 
points, and at most one 2-cycle. 


You saw in the main text that a fixed point a can be classified according to Recall that a fixed point a is 
the value of f’(a), the gradient of the graph of f at (a, f(a)). The gradient attracting if | f’(a)| < 1; 
formula for the function f(z) = x? +c is f’(x) = 22; in particular, for a repelling if | f'(a)| > 1; 

fixed point a, we have f’(a) = 2a. indifferent if | f’(a)| = 1. 


Similarly, a 2-cycle a,b of f can be classified according to the value of the A 2-cycle a,b is 


gradient product f’(a)f’(b); this is equal to the gradient of the graph of attracting if | f’(a)f’(b)| < 1; 
the composite function f o f at a, and also at b. For a 2-cycle a,b of the repelling if IF (a) FC )) > 1; 
function f(x) = x? +c, we have f’(a)f’(b) = (2a)(2b) = 4ab. indifferent if | f"(a) f'(b)| = 1. 


The first activity introduces the Mathcad worksheet that you will be using. 


Remember to create your own 
working copy of the Mathcad 
file. 


The variables M and p will be 
introduced in Activities 4.2 
and 4.5, respectively. You 
should not alter the values of 
these for now. 


Remember that to set 
variables in Mathcad we 
define them using the symbol 
:= , which can be obtained 
from the appropriate button 
on the ‘Calculator’ toolbar or 
by typing : (a colon, given by 
[Shift];). The symbol = 
evaluates variables in 
Mathcad. 
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Activity 4.1 Graphical iteration on the computer 


Open Mathcad file 221B1-01 Iterations of f(x)=x*2-+c, and move to 
page 2 of the worksheet. Here N iterations of the function f(x) = x7 +c 
are carried out, with initial term x9. The resulting sequence z,, is 
displayed on a graphical iteration diagram, and the last eleven terms 
calculated are listed in a table. Values for the fixed points, the 2-cycle, and 
the associated gradients and gradient product are displayed (to three 
decimal places) underneath the table. 


The page is set up with c = 0, so the function is initially f(x) = 2’. If you 
look at the fixed point and 2-cycle information, then you will see that 
Mathcad calculates that this function f has two fixed points, 0 and 1. The 
values of the gradients show that 0 is attracting whereas 1 is repelling. 
The 2-cycle variables indicate that f has no 2-cycle. 


The page is also set up with N = 10 and zo = 0.8. In parts (a) to (d) 
below you are asked to try other values for these two variables. 


(a) Describe the long-term behaviour of the iteration sequence with initial 
term rp = 0.8. Then set xo to each of the values —0.8, 0 and 1.1 in 
turn, and in each case describe the long-term behaviour of z,. 


(b) Set 2 9 = 0.999, which is a number just less than the repelling fixed 
point. Then set N = 100, and describe the long-term behaviour of z,. 


(c) With N still set to 100, set x) = 1.001, which is a number just greater 
than the repelling fixed point. Try to explain why Mathcad behaves as 
it does. 


(d) With zp still set to 1.001, set N = 15, and describe the long-term 
behaviour of z,,. 


Comment 


(a) When zo = 0.8, the sequence converges to the fixed point 0, with a 
staircase pattern of convergence. The same is true when 2p = —0.8. 


When zy = 0, the sequence converges to the fixed point 0. In fact, 
since 0 is a fixed point, each term of the sequence is equal to 0, which 
is why no iterations can be seen on the diagram. 


When zp = 1.1, the sequence tends to infinity. Only the first few 
iterations appear on the diagram, but you can see from the table that 
the terms of the sequence quickly become very large. 


When zp = 0.999, ten iterations are insufficient to indicate clearly the 
long-term behaviour of x,. When the number of iterations is increased 
to 100, the diagram shows the sequence approaching the fixed point 0. 
The table provides further evidence for this behaviour, since it shows 
that the terms Zoo tO Xj99 all have value 0 to three decimal places. 


When zp = 1.001 and N = 100, Mathcad marks f(z,) (above the 
graph) in red, and neither the table for x, nor the iteration diagram 
appears. Clicking on f(z,,) reveals an error message which shows that 
some of the terms of the sequence are too large for Mathcad to cope 
with. This suggests that the sequence is unbounded. To obtain a clear 
idea of the behaviour of the sequence, it is helpful to reduce the 
number of iterations, as suggested in part (d). 
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The terms in the table quickly become very large, which suggests that 
x, tends to infinity; this behaviour is also suggested by the diagram. 
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SECTION 4 Iterating real functions with the computer 


Mathcad notes 


© ‘The error that arose in part (c) occurs when the result of a calculation 
exceeds in magnitude the largest positive number that Mathcad can 
handle, which is about 102°’. 


© The expressions entered on the graph axes to plot the graphical 
iteration diagram have been ‘hidden’. This has been done so that the 
graph can fit on the page! The expressions on the axes of any graph 
can be hidden by clicking in the graph to select it, then choosing 
Graph p X-Y Plot... from the Format menu, selecting the ‘Traces’ 
tab and clicking in the check box to ‘Hide Arguments’. 


© The ‘staircases’ and ‘cobwebs’ are drawn on the graph by plotting 2,44 
against ©, using the trace type ‘step’ (see Figure 4.1). 


Mathcad plots the sequence of points (x,, 7,41), for k = 0,1,2,..., 
and joins each point (x,,%%41) to (p41, 2e42) using a horizontal line 
segment ending at (%,41, £441) followed by a vertical line segment. 

A separate trace is required to draw the first (vertical) line in the 
diagram from (9,0) to (a, #1). This is done by plotting x,-(k > 0) 
against ®p using the default trace type ‘lines’. Mathcad assigns the 
inequality (k > 0) in this expression the value 0 when it is false (for 
example, when k = 0) and 1 when it is true (for example, when k = 1). 


In the next activity you will begin to investigate how the sequences 
obtained by iterating the function f(x) = x? + c change as ¢ is varied. So 
that the investigation is conducted in a systematic manner, you should 
keep the initial term x9 set to the same value from now on. A suitable 
value for 2» is 0, and this is the value that will be used in the remainder of 
the investigation. Thus you will be investigating how the behaviour of the 
iteration sequence 


Re 23 =e +e iH =0,7,2,.:.) (4.1) 


varies as c is varied. Throughout Subsections 4.1 and 4.2, the notation zx, 
will always mean a sequence of the form defined in equation (4.1). 


For some values of c, the long-term behaviour of the sequence x,, may not 
be immediately obvious. As you work through the activity, you will be 
introduced to some ways in which you can adjust the values of variables on 
the Mathcad page to help you determine the long-term behaviour. 


Activity 4.2 Exploring iteration sequences 


In each of parts (a) to (e) you are asked to determine the long-term 
behaviour of the sequence z,, for a particular value of c, and then for two 
further values of c close to the first value. All of these values are marked 
on the number line in Figure 4.2, overleaf. 


An ‘F’ has been marked above the number 0 on the number line, to 
indicate that when c = 0 the sequence tends to a fixed point. As you work 
through the activity, record your findings for the other values of c in a 
similar manner. Use a ‘2’ if the sequence seems to tend to a 2-cycle, a ‘3’ 
for a 3-cycle, and so on. If the sequence seems to tend to infinity, then use 
the symbol ‘oo’. If it seems to be chaotic, that is, if it does not appear to 
settle down to any steady long-term behaviour, then use the letter ‘X’. 


First ensure that x 9 = 0. It is convenient to set the number of iterations 
and graph scale initially as follows: N = 10, s1 = —2 and s2 = 2. 


Remember that Mathcad 
notes are optional. 


Double-clicking in the graph 
gives an alternative approach. 


(2e-41, Tk+2) 


(Tk, 241) (Te+1, Se+1) 


Figure 4.1 Trace type ‘step’ 


You should still be working 
with Mathcad file 221B1-01. 


You looked at the case c = 0 
in Activity 4.1(a). 


Notice, however, that this 
sequence appears to be 
bounded: all terms seem to lie 
between about —1.8 and 1.5. 
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Figure 4.2 Values of c on the number line 


(a) Set c= 0.6. 


Record (on Figure 4.2) the long-term behaviour of the iteration 
sequence. Then repeat the process with c = 0.7 and c = 0.5. 


(b) Set c= 0.1. 
You may find it helpful to rescale the diagram, so that you can see the 
behaviour near the fixed point in more detail; try setting s1 = —0.2 
and s2 = 0.2. 
Record the long-term behaviour of the sequence. Then repeat the 
process with c = —0.1. 


(c) Reset s1 = —2 and s2 = 2. Set c= —0.6. 


You may find it helpful to increase the number of iterations and 
rescale the diagram; try N = 100, s1 = —0.8 and s2 = 0.2. 


Record the long-term behaviour of the sequence. Then repeat the 
process with c = —0.5 and c = —0.7. When c = —0.7, you may find it 
helpful to increase N still further; try N = 200. 


(d) Reset sl = —2, s2 = 2 and N = 10. Set c= —1.2. 


The long-term behaviour is not clear from the first ten iterations. 
However, the fixed-point and 2-cycle information below the graph 
shows that f has two repelling fixed points and an attracting 2-cycle, 
so the sequence may tend to this 2-cycle. To check this, increase the 
number of iterations; set N = 100. 


The values in the table do indeed appear to be settling down. To see 
the long-term behaviour clearly on the diagram, you can arrange for 
only the later iterations to be plotted. The graph is set up to omit the 
first M iterations from the plot, so you just need to increase the value 
of M; try M = 50. 
Record the long-term behaviour of the sequence. Then repeat the 
process with c = —1.1 and c= —1.3. 

(ce) Reset N = 10 and M =0. Set c= -1.8. 
Again, the long-term behaviour of the sequence is not clear. The 
fixed-point and 2-cycle information shows that although f has two 
fixed points and a 2-cycle, these are all repelling. 
Increase the value of N to see if the behaviour of the sequence seems 
to settle down. Try N = 100, to begin with, and then try increasing N 
still further, say to 200, and then to 300. You should find that new 
construction lines appear each time you increase NV, so it seems that 
the sequence does not settle down. It appears to be chaotic. 
Record the long-term behaviour of the sequence. Then repeat the 
process with c = —1.7 and c= —1.9. 


Solutions are given on page 34. 


SECTION 4 Iterating real functions with the computer 


In Activity 4.2 you were introduced to some ways in which you can adjust 


the values of variables in the Mathcad worksheet to help you find the 


long-term behaviour of an iteration sequence. These methods will be useful 


in the next activity. It is not possible to state an infallible strategy, but 
the list of suggestions in the box below should be a useful reference. 


Identifying the long-term behaviour of an iteration sequence 


The graphical iteration diagram, the table of terms, and the 
fixed-point and 2-cycle information can all help you to identify the 
long-term behaviour of an iteration sequence. Try the following 
suggestions. At any stage it may be helpful to rescale the diagram by 
adjusting the values of s1 and s2. 


To explore whether the sequence tends to a fixed point 


© Check that f has a fixed point, and that it is attracting (check 
that |f’(a)| <1, where a is the fixed point). 


© Increase the number JN of iterations until the sequence settles 
down and appears to tend to the fixed point. 


To explore whether the sequence tends to a 2-cycle 


© Check that f has a 2-cycle, and that it is attracting (check that 
|f'(a) f'(b)| < 1, where a and b are the members of the 2-cycle). 


© Increase the number of iterations until the sequence settles down. 
Then increase the number M of early iterations omitted from the 
plot until the remaining construction lines form a clear square. 


To explore whether the sequence tends to a p-cycle (p > 2) 


Increase the number of iterations until the sequence settles down. 
Then increase the number of early iterations omitted from the plot 
until only the more settled iterations are plotted. If possible, find the 
number of members of the cycle by counting the vertical construction 
lines in the diagram. 


To explore whether the sequence is unbounded 

Increase the number of iterations until the table shows that the terms 
become very large and positive or very large and negative, or the 
calculation for f(z,,) gives the error ‘Found a number with a 
magnitude greater than 10307 while trying to evaluate this 
expression. ’. 


To explore whether an iteration sequence is chaotic 

Increase the number of iterations (to, say, 100, then 200, then 300) 
and check that the sequence does not seem to settle down (check that 
new construction lines appear on the diagram each time the number 
of iterations is increased). 


From the results obtained in Activity 4.2, it appears that between the 
numbers —2 and 1 there may be a range of values of c for which the 
sequence x, tends to a fixed point, and a range of values of c for which it 
tends to a 2-cycle. In the next activity you are asked to find these ranges 
more precisely. 


The sequence has ‘settled 
down’ when increases to the 
number of iterations do not 
produce any visible new 
construction lines (assuming 
that no large area of the 
diagram is already completely 
covered by construction 
lines! ). 


If some members of the cycle 
are close together, then the 
graphical iteration diagram 
may not be clear enough to 
allow you to count them. 
Rescaling parts of the 
diagram may help. 


An apparently chaotic 
sequence may in fact tend to 
a cycle of very large order, 
but we ignore this possibility. 


You should still be working 
with Mathcad file 221B1-01. 


é ae 
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Figure 4.3 The case c= 0.5 


You should still be working 
with Mathcad file 221B1-01. 
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Activity 4.3 Attracting fixed points and attracting 2-cycles 


By testing further values of c in the Mathcad worksheet, try to determine 
the ranges of values of c between —2 and 1 for which the sequence x, 
behaves as follows: 


(a) x, tends to a fixed point; 
(b) x, tends to a 2-cycle. 


It may help to record the results of these tests with your earlier results in 
Figure 4.2. 


Solutions are given on page 34. 


In Activity 4.3 you saw a range of values of c for which the sequence x, 
tends to a fixed point, and a range of values of c for which it tends to a 
2-cycle. We can also state a range of values of c for which the sequence zr, 
tends to infinity. As c increases, the graph of y = 2* + c moves vertically 
upwards. For c > 0.25, all parts of the graph lie above the line y = z, and 
graphical iteration shows that the iteration sequence with initial term 

ro = 0 tends to infinity. This is illustrated by the graphical iteration 
diagram for c = 0.5 in Figure 4.3. Thus, for all c in the open interval 
(0.25, 00), the sequence z,, tends to infinity. 


In Activity 4.2, you saw that for c = —1.3 the sequence x, tends to a 
4-cycle, and you may have come across other values of c that give 4-cycles 
as you worked on Activity 4.3(b). If you wished, you could now go on to 
try to determine a range of values of c for which the sequence z,, tends to a 
4-cycle, and then perhaps go on to look at other values of c. However, for c 
less than about —1.4, the way in which the behaviour of the sequence 
varies as c varies is very complicated. Chaotic behaviour is possible, but 
attracting cycles can also occur, as you will see in the next activity. 


Activity 4.4 Attracting p-cycles 


Set sl = —2, s2 = 2, N = 100 and M = 50. 


(a) Set c to each of the following values in turn, and in each case state the 
long-term behaviour of the iteration sequence 2,: 


— 1,39) S475 2-257 HATE 


(b) Set c = —1.4, and try to decide whether the iteration sequence is 
chaotic or tends to a p-cycle. 


Solutions are given on page 34. 


Comment 


Once you have found a p-cycle that appears to be attracting, it is possible 
to check this classification of the p-cycle by calculating the corresponding 
gradient product; see the main text. 


SECTION 4 Iterating real functions with the computer 


The final, optional, activity in this subsection uses the Mathcad worksheet 

to illustrate the fact that if the function f has a p-cycle for some value 

of p, then the members of the p-cycle are all fixed points of the composite 

function f?. Recall that f? denotes the 
pth iterate of f, that is, the 


When the variable p in the Mathcad worksheet is set to an integer greater 
composite function 


than 1 (and not too large), the graph of f? is added to the graphical fo fio---o f, obtained when 
iteration diagram. f is applied p times; see the 
main text. 


Activity 4.5 p-cycles and the function f” (Optional) 


pet at = —2, s2 = 2, N = 100 and 1 = 50. You should still be working 
(a) Set c= —1.2 and p= 2. Observe from the graph of f? that the with Mathcad file 221B1-01. 
members of the 2-cycle are attracting fixed points of f?. 


(b) Repeat the process for the following values of ¢ and p: You have already considered 


these values of c, in 
tS -iap=s ¢= 145, p=6, cae p= 3. Activities 4.2 and 4.4. 


Mathcad notes 


The function f? is displayed on the graph only when p = 2,3,4,.... Its 
suppression when p is 0 or 1 is achieved by multiplying the expressions 
entered on the graph axes to plot the graph of f? by (p > 1). Mathcad 
assigns the expression (p > 1) the value 0 when it is false and 1 when it is 
true. 


Now close Mathcad file 221B1-01. 


You may have wondered how we obtained the values of c given in 
Activity 4.4. It is not easy to find such attracting cycles by testing 
‘random’ values of c; if you experimented using the Mathcad worksheet 
with values of c less than about —1.4, then you probably found that you 
usually obtained a sequence whose behaviour appeared to be chaotic. The 
next subsection shows how to find values of c that give cycles. 


4.2 The bifurcation diagram 


The diagram in Figure 4.4, overleaf, provides an overview of the changes in 

the long-term behaviour of the sequence x, as c is varied. The diagram was Recall that the sequence z,, is 
generated, using Mathcad, by considering 1501 values of c, equally-spaced defined by the recurrence 

on the horizontal axis. For each of these values of c, the iteration sequence System in equation (4.1). 


Pep P GES; Sj DSe 
was calculated. The first 200 terms of this sequence were then ignored, and 
the remaining terms displayed on the graph, by plotting each of the points 

(c, L200), (c, £201), a ie £300); 
as a small dot. 


If you pick a value of c, and imagine a vertical line drawn through the 
point on the horizontal axis of the diagram corresponding to that value 
of c, then the points where that vertical line intersects the graph give an 
indication of the long-term behaviour of the sequence z,,. 
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A single intersection point 
occurs when the plotted 
terms of the sequence are 
very close together. ‘Two 
intersection points occur 
when each plotted term is 
close to one of two values. 
Three or more intersection 


points occur in a similar way. 


To ‘bifurcate’ is to fork into 
two. 
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If there is a single point of intersection, then this suggests that the 
sequence tends to a fixed point. Two intersection points suggest that it 
tends to a 2-cycle, 3 points suggest a 3-cycle, and so on. Many points 
suggest chaotic behaviour. You can see from the diagram, for example, 
that values of c in the interval (—1.25, —0.75) seem to yield sequences 2, 
that tend to a 2-cycle. This agrees with the solution to Activity 4.3(b). 


2 
“9 -cat8scrbbie eae la poetic 08 oc Dee ee 
Figure 4.4 A bifurcation diagram for functions of the form f(x) = r* +c 


The diagram illustrates that there are certain key values of c where the 
behaviour of the sequence z,, changes. For example, as c decreases through 
—0.75, the sequence stops tending to a fixed point, and begins to tend to a 
2-cycle. Similarly, as c decreases through —1.25, the sequence stops 
tending to a 2-cycle, and begins to tend to a 4-cycle. At such values as. 
each member of the cycle ‘forks’ into two cycle members, and we say that 
a bifurcation occurs. The graph in Figure 4.4 is known as a bifurcation 
diagram. 


At each bifurcation the number of members of the cycle doubles, since 
every member of the cycle ‘forks’ into two. For example, as c decreases 
from the value 0, sequences tending to fixed points give way to sequences 
tending to 2-cycles, which give way to sequences tending to 4-cycles, and 
so on. Similarly, as c decreases from the value —1.75, sequences tending to 
3-cycles give way to sequences tending to 6-cycles, and so on. There are in 
fact infinitely many ‘windows’ in which the sequences tend to cycles and 
doubling of the number of members of the cycles occurs; in each of them 
the bifurcation points become progressively closer together as c decreases. 


Outside these windows, sequences show bounded chaotic behaviour. For 
example, when c = —1.8, the sequence appears to be chaotic, with all 
terms lying between approximately —1.8 and 1.5; this confirms what you 
saw in Activity 4.2(e). You can see from the bifurcation diagram that as c 
decreases, the interval of values within which terms of the sequence lie 
gradually expands. 


In the next activity you will use Mathcad to explore the bifurcation 
diagram. You will be asked to try to find a value of c that gives an 
iteration sequence which tends to a 5-cycle. 


SECTION 4 _ Iterating real functions with the computer 


Activity 4.6 The bifurcation diagram. 


Open Mathcad file 221B1-02 Overview of iterations of f(x)=x’2-+c. 
This worksheet draws a bifurcation diagram for functions of the form 

f(x) = x? +c, with initial term x9 = 0, for values of c from C1 to C2. You 
may have to wait for some time (perhaps 30 seconds) while Mathcad 
performs the calculations needed to draw the graph. 


The bifurcation diagram is produced by taking V + 1 equally-spaced values 
of c. For each of these values of c, N iterations are carried out, and the 
first M of these are omitted from the plot. The worksheet is set up with 
=a =. 300 and At = 200. 


You can see parts of the diagram in more detail by altering the values of 
the horizontal axis limits Cl and C2. You can also increase the value of V; 
this produces a clearer diagram, but the calculations take longer. 


Try to find a value of c with an attracting 5-cycle, by proceeding as follows. 


(a) First set C2 = —1.5, so that you can see the part of the diagram 
corresponding to values of c between —2 and —1.5 in more detail. 


(b) Look at the part of the diagram between —1.65 and —1.6. You should 
see a narrow window that appears to include values of c for which the 
sequence z, tends to a 5-cycle. Set Cl = —1.65 and C2 = —1.6 to 
confirm this observation, and choose a value of c that gives a 5-cycle. 


(c) If you wish, test your value of c in Mathcad file 221B1-01. 
A solution to part (b) is given on page 34. 


Comment 


© You can interrupt a Mathcad calculation by pressing [Esc], the 
escape key, and then clicking on OK in the resulting option box. 
Mathcad will recalculate when you next make a change. 


© By default, Mathcad operates in ‘automatic calculation mode’, but this 
can be inconvenient where more than one input change is to be made 
before a time-consuming recalculation is required. In order to switch to 
manual calculation mode, which disables automatic calculation, select 
Automatic Calculation from the Math menu. (When you have 
done this, the tick mark beside Automatic Calculation in the menu 
disappears, and the word ‘AUTO’ disappears from the status bar in 
the bottom right corner of the Mathcad window.) Once in manual 
mode, you can calculate results when you choose, either by selecting 
Calculate from the Math menu, or by pressing the [F9] function key. 


Mathcad notes 


© The graph is produced by plotting x; against c; using the trace type 
‘points’. The subscripted variable x;,,, is entered in the usual way. 
Either use the ‘x,’ button on the ‘Matrix’ toolbar or type [ (left 
square bracket), then separate the subscripts i and n with a comma. 
Thus you could obtain z;,, on the screen by typing x[i,n. 


© A Mathcad graph can display at least 150000 individual points. If you 
try to plot a graph with more points than this, then an error may 
occur. (No graph is drawn and the graph box is highlighted in red — 
clicking on it reveals the error message ‘Can’t plot this many points. ’.) 


Now close Mathcad file 221B1-02. 


The bifurcation diagram in 
Figure 4.4 was produced by 
taking V = 1500, N = 300 
and M = 200. 


If you seek to change both C1 
and C’2, then you will find 
that Mathcad starts to 
recalculate after the first 
change has been made. The 
Comment below gives ways to 
deal with this problem. 


More information about this 
technique, and details of how 
to interrupt and resume 
calculations, are provided in 
A Guide to Mathcad. 


In order to return to 
automatic mode, follow the 
same procedure, whereupon 
the tick mark reappears. 


Lin is the notation used in 
the worksheet for the nth 
term of the sequence obtained 
using the 7th value of c, 
denoted by c;. 
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If you have studied 

Chapter B1 of MST121, then 
you may recognise that this 
sequence is associated with 
the logistic recurrence 
relation, and you will have 
seen the diagram in Figure 4.5 
in connection with this. 


Essentially the same 
bifurcation diagram occurs 
for various families of 
iteration sequences, generated 
both by quadratic functions 
and by other functions. 


You can check this, if you 
wish, by expressing z,, in 
terms of x, and substituting 
for to, Z, and 7,.1.m 
equation (4.2). 
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4.3 Other quadratic iteration sequences (Optional) 


Other families of quadratic functions also show complicated behaviour 
when they are iterated. Consider, for example, the family of iteration 
sequences 


ds = 1S, tas = SF Peal — ty) ee) (4.2) 


obtained by iterating the quadratic function f(x) = «+ ra(1-— <2), where r 
is a parameter. Different values of r determine different iteration sequences 
from this family, in the same way that different values of c determine 
different iteration sequences from the family in equation (4.1). Figure 4.5 
shows a bifurcation diagram for this new family, with r between 1 and 3. 
The structure of the graph appears essentially the same as that of the 
bifurcation diagram in Figure 4.4, although it is reversed and distorted. 


Lie 


0.5 


0 


l la 46 ee 2 p= eee) See >: ea 3 


Figure 4.5 A bifurcation diagram for functions of the form f(x) =x +rz(1— 2) 


The similarity of the bifurcation diagrams for these two families can be 
explained as follows. Suppose that the sequence z,, satisfies equation (4.2). 
Consider the new sequence 


x! = —ran + 4(1+r) (n=0,1,2,...). 


This sequence z’, is related to x, by scaling (by the factor —r) and shifting 
(by adding $(1+1r)). Thus both the sequences z, and x, have the same 
type of long-term behaviour (for example, if x, tends to a 3-cycle, then 2’, 
tends to a 3-cycle). Now, it can be shown that the new sequence 7, 
satisfies 


, = US, iy te Pe = Cee, «3, (4.3) 


where c = +(1—r*). Equation (4.3) is of the same form as equation (4.1), 
except that the initial term has changed from 0 to 0.5. It can be shown 
that this change of initial term does not affect the long-term behaviour of 
these sequences. Thus the long-term behaviour seen on a vertical slice 
through the graph in Figure 4.5, for some particular value of r, is of the 
same type as that seen on the vertical slice through the graph in Figure 4.4 
for the corresponding value c = 4(1— Tr’). 


Now as the parameter r takes values increasing from 1 to 3, the parameter 
c = (1 —r?) takes values decreasing from 0 to —2. Thus the patterns in 
Figure 4.5 occur in the reverse order to those in Figure 4.4. They are also 
more ‘squashed together’ at one end because 4(1 — r*) decreases 
progressively faster as r increases from 1 to 3. Other distortions occur 
because the scaling and shifting factors relating x, to x}, vary as r varies. 


SECTION 4 Iterating real functions with the computer 


4.4 Sequences obtained by iterating real functions 


In this subsection you will use Mathcad file 221B1-03, which allows you to 
explore the sequences obtained by iterating any real function. 


If we wish to study the long-term behaviour of the sequences obtained by 
iterating a given real function f, then it is useful to find and classify any 
fixed points and 2-cycles of f. The Mathcad worksheet uses the Mathcad 
solve block to find approximate solutions to the fixed point equation and 
the 2-cycle equation. The solve block can be used whether or not these 
equations have solutions given by formulas. You have to provide Mathcad 
with a ‘guess’ for a solution; it then uses an iterative method to calculate a 
sequence of values that become progressively closer to a solution. Mathcad 
usually obtains a solution accurate to several decimal places; if it cannot 
find one, then it registers an error. Different initial guesses may yield 
different solutions. 


To classify the fixed points and 2-cycles of a function f, we have to be able. 


to find the gradient of the graph of f at the fixed points, and at members 
of the 2-cycles. The Mathcad worksheet uses a built-in feature of Mathcad, 
the d/dz operator, to find an approximate value for the gradient, here 
denoted by the expression Df(x). You will learn much more about this 
topic in Block C; for now, just accept the values that Mathcad provides! 


In the next activity you are asked to explore iterations of a particular real 
function, using the following strategy. 


Exploring iteration sequences using Mathcad file 221B1-03 
Edit the definition of the function f, on page 2, as required. 


To find and classify all the fixed points of f (page 2) 


1. If necessary, rescale the graph to show all the fixed points of f. 
(Adjust the values of the axis limits s1 and s2.) 


2. Use the graph to estimate the approximate value of a fixed point 
of f, and set the ‘guess’ value x to this value. Read off the more 
accurate value a given for the fixed point. Use the value of the 
gradient f’(a) to classify the fixed point. 


3. Repeat the instructions in step 2 as many times as is necessary to 
find and classify all the fixed points of f. 


To find and classify all the 2-cycles of f (page 3) 

1. If necessary, rescale the graph to show all the fixed points of the 
composite function f o f. Identify the fixed points of f o f that 
are not fixed points of f (these are the points where the graph of 
fo f meets the line y = x but the graph of f does not). These 
points pair off into 2-cycles of f. 


2. Use the graph to estimate the approximate value of a fixed point 
of f o f that is not a fixed point of f, and set the ‘guess’ value x 
to this value. Read off the more accurate value a, and the 
corresponding value b = f(a), for the 2-cycle a,b of f. Use the 
value of the gradient product f’(a)f’(b) to classify the 2-cycle. 


3. Repeat the instructions in step 2 as many times as is necessary to 
find and classify all the 2-cycles of f. 


To identify the long-term behaviour of a sequence (page 4) 


Use the suggestions given in the box on page 9 to help pa to find the 
long-term behaviour of the iteration sequence. 


A strategy for using this file 
is given in the box below. 


If you have done the 
computer work for 

Chapter A3 of MST121, then 
you will have met the solve 
block. It is also covered in 

A Guide to Mathcad. 


In Block C, you will see that 
the gradient f’() of the 
graph of a function f at 

(x, f(x)) is also represented 
by the notation 


d 
ay \f(e)): 


called the derivative of f at z. 


This strategy should be read 
while you are doing 


Activity 4.7 overleaf. 


The value a is indicated on 
the graph by a small black 
box at the point (a,a). If a is 
not the value that you wished 
to find, then try to set the 
‘guess’ value x more 
accurately. 


The values a and 6 are 
indicated on the graph by 
small black boxes at the 
points (a,a) and (6, b), joined 
by dashed black lines. If a is 
not the value that you wished 
to find, then try to set the 
‘guess’ value x more 
accurately. 


The variables NV, M, s1, s2 
and p have the same roles 
here as in Subsection 4.1. 
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The prime symbol is obtained 
in Mathcad by pressing the 
left quote key. 
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Activity 4.7 Iterating real functions 
Open Mathcad file 221B1-03 Iterating real functions. The function 


Ar 
Re SZ (1 + 2} 
is already entered for you on page 2 of the worksheet. Use the worksheet 
and the strategy on page 15 to help you to do the following. 


(a) Find all the fixed points of f, and classify them as attracting, repelling 
or indifferent. 


(b) Find all the 2-cycles of f, and classify them as attracting, repelling or 
indifferent. 


(c) For each of the following initial terms xo, find the long-term behaviour 
of the sequence obtained by iterating f with that initial term: 


5 = Sa, Sg 
Solutions are given on page 34. 


Mathcad notes 


© Both the solve block (used to find a solution of an equation) and the 
d/dx operator (used to calculate a gradient) use numerical methods to 
obtain approximations to the exact solution and gradient, respectively. 
The error messages ‘No solution was found. ...’ (for the solve block) 
and ‘Can’t converge to a solution.’ (for the d/dzx operator) indicate 
that the numerical method has failed. In the case of the solve block it 
is worth trying different initial guesses, though there may in fact be no 
solution. In the case of the d/dz operator, you could try to find the 
gradient at a nearby point. 


© 


The accuracy of the solve block is controlled by the built-in variable 
TOL (Math menu, Options..., ‘Built-In Variables’ tab, Convergence 
Tolerance (TOL)). Mathcad looks for a solution until the difference 
between two successive approximations is less than or equal to TOL. 
By default, TOL = 0.001, but for this file it is set to 0.000001. For the 
d/da operator, the value given is usually accurate to 7 or 8 significant 
figures, irrespective of the value of TOL. 


© The choice of the name Df for the gradient function of f in the 
Mathcad worksheet is a pragmatic one. The usual notation, f’ 
(‘f prime’ or ‘f dash’), could have been used, but the prime symbol is 
hard to see when placed to the right of an ‘f’ in Mathcad. 


Now close Mathcad file 221B1-08. 


Chapter B2, Section 5 
Visualising affine transformations 


Mathcad provides a convenient means by which you can see the effects of 
linear and affine transformations on the unit square and unit grid. 


5.1 Linear transformations 


In the first activity you will explore linear transformations of the plane. 
You will use a Mathcad worksheet that allows you to enter any 2 x 2 
matrix A and see the effect of the linear transformation f(x) = Ax on the 
unit square and unit grid. 


You can take f to be one of the basic linear transformations by setting A 
equal to one of the following matrices, expressed in ‘Mathcad notation’. 


Mathcad matrix notation for basic linear transformations This notation is used in all 
the Mathcad worksheets for 


Linear transformation Mathcad matrix notation Chapters B2 and B3. 


Rotation about the origin 
through angle 0 


Reflection in the line through 
the origin that makes angle 6 
with the positive x-axis 


Scaling with factor a in the 
x-direction and factor 6b in the 
y-direction 


Uniform scaling with factor a 


xz-shear with factor a 


y-shear with factor a 


For example, the matrix of a rotation about the origin through angle 7/2 
anticlockwise can be specified in the worksheet by entering ‘R(Z)’. Note 
that our Mathcad notation for rotation and reflection matrices differs 
slightly from the notation used in the main text. In the Mathcad notation, 
parameters appear in brackets (for example, R(5)), whereas in the main 
text they appear as subscripts (for example, R, 2). The Mathcad notation 
is used not only in the Mathcad worksheets but also in this computer book. 


In the first worksheet, the side of the unit square that lies along the x-axis 
is marked with a filled-in triangle, as shown in Figure 5.1(a), overleaf. The 
image of this triangle is plotted on the image square, so that you can tell 
which vertices of the square have been mapped to which vertices of the 
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See Chapter B2, Section 2. 


These matrix techniques are 
used in the computer work for 
Chapter B2 of MST121. 


The determinant of A is 
denoted by |A| in Mathcad. 


The edges of the graph box 
may be obscured in places by 
the unit grid. 


For parts (g) and (h), define 
A :=M. The matrix in 

part (g) is already entered for 
you; edit M to obtain the 
matrices in part (h). 
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image. This lets you see whether orientation has been preserved or 
reversed. For example, Figure 5.1(b) shows the image of the unit square 
under a rotation through 7/2 anticlockwise about the origin, and 
Figure 5.1(c) shows its image under reflection in the y-axis. 


(a) (b) (c) : 


Figure 5.1 The unit square and two images 


You saw in the main text that you can predict whether a linear 
transformation f preserves or reverses orientation by considering the 
determinant of its matrix A. The transformation f preserves orientation if 
det A > 0 and reverses orientation if det A < 0. If det A = 0, then f isa 
flattening, which destroys orientation. The determinant of A also tells you 
the effect of f on areas: these are scaled by the factor | det Al. 


Nene een 


Activity 5.1 Images of the unit square and unit grid 


Open Mathcad file 221B2-01 Linear transformations. Page 2 of the 
worksheet describes the basic techniques for creating, editing and 
calculating with matrices in Mathcad. Page 3 sets up and explains the 
notation used in the worksheet. 


Move to page 4. Here a 2 x 2 matrix A is defined, and its determinant 
calculated. Graphs display the unit square and unit grid, and their images 
under the linear transformation f(x) = Ax. You can set A to be a matrix 
representing one of the basic linear transformations, using the notation 
explained before this activity. Alternatively, you can set A to be the 
general 2 x 2 matrix M. The page is set up with A = R (2). 


The graphs can be rescaled by changing the value of s, which is defined 
near the top of the page; both axes of both graphs are scaled from —s to s. 


Set A to be each of the following matrices in turn. In each case, use the 
value of the determinant of A to predict whether the transformation 
preserves, reverses or destroys orientation, and whether it decreases, 
preserves or increases area, and try to confirm your predictions by looking 
carefully at the effect of the transformation on the unit square. 


a) Rotations: R(4), R(—7). 
b) Reflections: Q(4), Q(=). 
c). Scalings: S(2,1), $(0:5, 1), S(—1,1),.5(=1, <3}. 
Uniform scalings: U(2), U(0.5), U(—1), U(—0.5). 
e) x-shears: X(2), X(1), X(0), X(—1), X(—2). 

) y-shears: Y(2), Y(—1). 


4 1-6 
eoel nat 


ae 4 
Non-basic linear transformations: ( 9 


g) A flattening: ( 


SECTION 5 _ Visualising affine transformations 


Comment 


© Parts (a) to (g) illustrate the effects of basic linear transformations on 
orientation and area, which can be summarised as follows. 


Orientation-preserving transformations: rotations, shears, scalings Remember that a uniform 


with factors a, b of the same sign. scaling U(a) is a scaling, with 


: : : ; both factors equal to a. 
Orientation-reversing transformations: reflections, scalings with factors 


a, b of opposite signs. 


Area-preserving transformations: rotations, reflections, shears, scalings 
with factors a, b where |ab| = 1. 


Area-decreasing transformations: scalings with factors a, b where 
lab] < 1, flattenings. 


Area-increasing transformations: scalings with factors a, b where 
jab] > 1. 
Flattenings destroy orientation and decrease area to zero. 


© The linear transformation represented by the first matrix in part (h) 
preserves orientation and increases area; the second reverses 
orientation and decreases area. 


We now consider general linear transformations. For example, Figure 5.2 
shows the effect on the unit square and unit grid of the linear 


transformation m represented by the matrix M = ( : Zs ) 


Figure 5.2. Image under the linear transformation m 


This linear transformation is in fact the composite of a uniform scaling 
with factor 2, followed by an x-shear with factor 1, followed by a rotation 
through 7/2, as we now check using matrix multiplication. 


Recall that if f and g are linear transformations represented by matrices A See Chapter B2, Section 3. 
and B respectively, then the composite linear transformation go f is Remember that go f means 
represented by the product matrix BA. It follows that if f, g and h are ‘first apply f, then apply gq’. 
linear transformations of the plane represented by matrices A, B and C 

respectively, then the composite linear transformation ho go f is 

represented by the matrix CBA. 


If we take A to be the matrix representing a uniform scaling with factor 2, 
B to be the matrix representing an x-shear with factor 1, and C to be the 
matrix representing a rotation through 7/2, then 


SBA = (aod ea caditee. a1) Taye 09) 


as required. 


19 


The uniform scaling with 
factor 1, the x-shear with 
factor 0 and the rotation 
through angle 0 are all equal 
to the identity 
transformation. 


(a) (b) (c) 


CHAPTER B2 


Any linear transformation m of the plane can be expressed as a composite 
of at most three basic linear transformations. This shows that every linear 
transformation of the plane has a fairly simple geometrical interpretation. 
In particular, if m is not a flattening, then it can be expressed as the 
composite of a scaling, followed by an a-shear, followed by a rotation. 


Some of the basic linear transformations in this ‘decomposition’ of m may 
be the identity transformation, in which case m is a composite of two basic 
linear transformations, or is itself a basic linear transformation. 


To express a given linear transformation that is not a flattening as a 
composite of a scaling, followed by an x-shear, followed by a rotation, we 
consider its effect on the unit square. Figure 5.3(a) shows the unit square 
and Figure 5.3(d) shows its image under a typical linear transformation m. 
The image is a parallelogram with one vertex at the origin, which we call 
the target parallelogram. We use the word base to describe the side of the 
unit square that lies along the z-axis and the image of this side on the 
target parallelogram. These base sides are marked with filled-in triangles. 
Figure 5.3(b) and (c) illustrate stages in finding the decomposition, with 
the corresponding base sides marked similarly. 


Figure 5.3 Successive images of the unit square 


This rectangle has the same 
orientation and area as the 
target parallelogram. 


Here d is positive if the point 
moves to the right and 
negative if it moves to the 
left. 
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The process of finding the basic linear transformations that constitute the 
decomposition of m is described below. 


Decomposition of a linear transformation ™m 


Suppose that the target parallelogram under m has a base of length 1, 
and height h, as shown in Figure 5.3(d). The three basic linear 
transformations for the decomposition can be chosen as follows. 


1. First choose the scaling with matrix S(+l,h), where the + sign is 
used if m preserves orientation and the — sign otherwise. This 
maps the unit square to a rectangle with base of length / and 
height h; see Figure 5.3(b). 

Next take the x-shear that maps the rectangle onto a parallogram 
that can be rotated onto the target parallelogram; see 
Figure 5.3(c). 


Finally, take the rotation about the origin that maps the 
parallelogram onto the target parallelogram in Figure 5.3(d). 


In step 2, you can find the x-shear required by determining the signed 
distance, d say, that the x-shear moves a point on the side of the rectangle 
opposite the base. The shear factor is then d/h. 


In the next activity you will use a page of the Mathcad worksheet that is 
designed to help you carry out this decomposition process for a given 
linear transformation. Parts (a) and (b) of the activity provide examples, 
and in part (c) you are asked to work through the process yourself. 


SECTION 5 Visualising affine transformations 


Activity 5.2 Composite linear transformations 


Move to page 5 of the worksheet. Near the top of the page a matrix Mis —__ You should still be working 
defined, and the value of its determinant |M| is displayed. Graphs display | with Mathcad file 221B2-01. 
the unit square and unit grid, and their images under the linear 

transformation m(x) = Mx. The bases of the unit square and the target 

parallelogram (shown in blue) are marked with filled-in triangles. 


Beneath this, the page is set up so that you can define the matrix A of a 
scaling, the matrix B of an x-shear, and the matrix C of a rotation. The 
determinants of A, B, C and CBA are all displayed. Initially, we have 


A=S01), Bare, == ane. (5.1) Note that each of the 
matrices S(1,1), X(0) and 


The linear transformations represented by A, B and C are defined as f, g BO cectesents he dane 


and h, respectively. The images of the unit square and unit grid under f, 


transformation. 
go f,andhogof are displayed in successive graphs. Each of these graphs 
also shows the target parallelogram, in dashed blue lines. All the graphs 
can be rescaled by changing the value of s, defined near the top of the page. 
(a) The page is set up with 
0 -2 | = 
M = . 5): The decomposition of this 
matrix M was discussed 
In this case |M| is positive, so m preserves orientation. earlier in the subsection. 


Set each of the matrices A, B and C as shown below, in turn, and as 
you do so check that they achieve steps 1, 2 and 3, respectively, of the 
decomposition process on page 20. 


A =S(2,2), B=X(1), C=R(2). 


(b) Now reset each of the matrices A, B and C as in equations (5.1). 
Then edit the entries of the matrix M to give 


—j] -I 
Meo 4) 
In this case |M| is negative, so m reverses orientation. 


Set each of the matrices A, B and C as shown below, in turn, and as 
you do so check that they achieve steps 1, 2 and 3, respectively, of the 
decomposition process on page 20. 


A=S(-1,2), B=X(-1), C=R(0). 


(c) In each of parts (i) to (iv) below, first reset each of the matrices A, B 
and C as in equations (5.1), and set M to the given matrix. Then 
follow the decomposition process to express m as a composite of basic — For each of these matrices, 
linear transformations. Hence express each given matrix as a product — the parameters of the 
of matrices of basic linear transformations. required scaling and shear are 
2 9 = oe 7 oe 0 3 small integers or simple 
(i) ( ay ) (ii) ( eee | ) (iii) ( — ) (iv) ( =: fractions, and the parameter 
of the required rotation is an 
Solutions to part (c) are given on page 35. integer multiple of 7/2. 


Now close Mathcad file 221B2-01. 


ra 


In each of the given examples, 
A is the matrix of a basic 
linear transformation, and the 
vector a has integer 
components. 


However, it is not necessary 
to view this area in order to 
carry out the activities based 
on the worksheet, nor to 
understand the principles 
behind them. 
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5.2 Affine transformations 


In the next activity you will use Mathcad to visualise the effect of some 
affine transformations of the plane. You will be given the images of the 
unit square and unit grid under various ‘mystery’ affine transformations, 
and your task is to find the affine transformations that produce these 
images. 


Activity 5.3 Exploring affine transformations 


Open Mathcad file 221B2-02 Affine transformations. Page 2 of the 
worksheet defines the matrices representing the basic linear 
transformations, using the Mathcad notation. 


Move to page 3. The unit square and unit grid, and their images under a 
‘mystery’ affine transformation, are displayed near the top of the page. 
The page can be set up with any one of four different mystery 
transformations, and you can choose which one of these is applied by 
setting the variable T to 1, 2, 3 or 4. All of the graphs can be rescaled by 
changing the value of s, defined near the top of the page. 


You can define your own affine transformation f(x) = Ax + a by editing 
the definitions of the matrix A and vector a near the middle of the page. 
The unit square and unit grid, and their images under f, are displayed 
immediately below these definitions. The graph that displays the images 
under f also shows the image of the unit square under the mystery affine 
transformation, in dashed blue lines; this is the target. 


Set T to each of 1, 2, 3 and 4 in turn. In each case, try to identify the 
mystery affine transformation, and verify your answer by setting A and a 
and checking that you hit the target. 


Solutions are given on page 35. 


Comment 


You may wish to use the Mathcad page to see the effects of other affine 
transformations on the unit square and unit grid. If you set T’ = 0, then no 
target will appear on the graphs. 


Mathcad notes 


The definitions of the four mystery affine transformations are hidden in the 
columns of a matrix off to the right of page 3 of the worksheet. The area 
beyond the right-hand margin of a Mathcad page (which is marked by a 
solid vertical line) can be used just like the rest of the worksheet. It is 
divided into further pages, where you can place mathematical expressions, 
text, graphs and pictures. Some of the other Mathcad worksheets in the 
course also have material off the page to the right. 


You can view the pages in this area by using the horizontal scroll bar to 
move to the right. When printing a wide worksheet, you can choose 
whether or not the content to the right of the left-hand pages is to be 
printed. To do this, select Page Setup... from the File menu. Then, in 
the ‘Page Setup’ option box, tick or untick the check box for ‘Print single 
page width’, as required. If there is a tick here, then just the left-hand 
pages of the worksheet will be printed when Print... is selected 
subsequently; this is the default for MS221 worksheets. Otherwise, the 
whole width of the worksheet will be printed. 


SECTION 5 _ Visualising affine transformations 


Given any three points, in any order, there is a unique affine 
transformation f that maps the points (0,0), (1,0), (0,1) to these points, 
respectively. ‘The main text gave a method for finding f. The next, 
optional, activity allows you to explore the effect of affine transformations 
on the triangle with vertices (0,0), (1,0), (0,1), which we call the unit 
triangle. ‘This lets you check visually, using Mathcad, that an affine 
transformation found using the given method does indeed map the three 
points (0,0), (1,0), (0,1) to the required image points. 


Activity 5.4 Applying affine transformations to the unit 
triangle (Optional) 


Move to page 4 of the worksheet. Here you can define an affine 
transformation f(x) = Ax +a by setting the matrix A and vector a as 
described in the worksheet. Graphs display the unit grid and unit triangle, 
and their images under f. The coordinates of the vertices of the image 
triangle, and its area, are shown at the bottom of the page. 


Use the Mathcad page to verify the correctness of the affine 
transformations found in Example 4.1 and Activity 4.2 in the main text. 


Mathcad notes 


The ‘|x|’ button (available on both the ‘Calculator’ and ‘Matrix’ toolbars), 
whose keyboard alternative is [Shift]\ (shift and backslash), performs 
several roles in Mathcad. It gives the determinant of a matrix, the 
modulus (absolute value) of a real number, the magnitude of a vector, and 
the modulus of a complex number. 


The modulus of the determinant of a matrix A (written as | det A| in the 
main text) can be obtained by two applications of ‘|x|’, and appears 
as ||A||. 


Now close Mathcad file 221B2-02. 


See Chapter B2, Section 4. 


You should still be working 
with Mathcad file 221B2-02. 


You will study complex 
numbers later in the course. 


ae 


Chapter B3, Section 5 
Iterating linear transformations with the computer 


This notation is given on 
page 17. 


See Chapter B3, Section 4. 


Property (c) states that if ky 
is the ‘dominant eigenvalue’, 
then (2,, Yn) tends in the 
direction of the ‘dominant 
eigenline’ y = mz. 
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In this section you will use Mathcad to explore sequences of points in the 
plane obtained by iterating linear transformations. You will also see a few 
examples of sequences obtained by iterating affine transformations. The 
section ends with an optional subsection in which you can see some 
visually interesting plots obtained by an iteration process involving more 
than one affine transformation. 


The Mathcad notation for matrices of basic linear transformations used in 
the computer work for Chapter B2 will be used again in this section, both 
in the text and in the Mathcad files. 


5.1 Iterating linear transformations 


This subsection is concerned with sequences of points in the plane 
obtained by iterating linear transformations. In the main text, you saw 
some examples of such sequences where the matrix representing the linear 
transformation has two distinct non-zero eigenvalues, the case of a 
so-called generalised scaling. These have the following properties. 


Iteration properties of generalised scalings 

Let the linear transformation f be represented by a 2 x 2 matrix A 

that has two distinct non-zero eigenvalues k, and k2, with 

corresponding eigenlines @; and @). Let (xo, yo) be a point of R? and 

let (2, Yn) be an iteration sequence generated by A, with initial 

point (Xo, Yo): 

(a) (i) If k, > 0, then (x, yn) all lie on the same side of £2 as (Zo, Yo). 
(ii) If ky < 0, then (z,,y,) alternate between opposite sides of 05. 


(b) (i) If max{|ky|,|ko|} > 1, then the sequence moves away 
from (0,0). 
(ii) If max{|k,|, |k2|} < 1, then the sequence moves towards (0,0). 


(c) If |k,| > |k2| and (xo, yo) does not lie on an eigenline, then 


where m is the gradient of @;. 


The first activity involves sequences obtained by iterating non-uniform 
scalings. These are the ‘simplest’ generalised scalings. The scaling with 
factors a and b, where a # b, has eigenvalues a and b and the corresponding 
eigenlines are y = 0 (the z-axis) and x = 0 (the y-axis), respectively. 


For example, Figure 5.1 shows the first few points of the iteration sequence 
generated by A = $(1.5,1) and the initial point (2,1). It illustrates 
behaviour that can be predicted using the iteration properties of 
generalised scalings given in the box above. Since both eigenvalues are 
positive, all points of the sequence lie on the same side of the eigenline 

y = 0, and on the same side of the eigenline x = 0, as the initial point. 
Since max{|1.5|,|1|} > 1, the sequence moves away from (0,0). 


SECTION 5 Iterating linear transformations with the computer 


Since |1.5| > |1|, the dominant eigenvalue is 1.5, the dominant eigenline is 
y = 0, and the sequence tends in the direction of this line. 


It follows from property (b) of 
generalised scalings (see 
Chapter B3, Section 3) that 
this sequence remains the 
same distance from the line 

y = 0 and moves away from 
the line x = 0. 


Figure 5.1 Iteration sequence generated by a non-uniform scaling 


In Activity 5.1 you will see sequences obtained by iterating various 
non-uniform scalings. 


Activity 5.1 Non-uniform scalings 


Open Mathcad file 221B3-01 Iterating linear transformations. 
Page 2 of the worksheet defines the matrices representing basic linear 
transformations, using the Mathcad notation. 


Move to page 3. Here a 2 x 2 matrix A and an initial point (29, yo) are 


defined; the initial point is defined as the vector xp = ( S ). In this computer book, 
0 

A graph displays the first point, as a black box symbol, and N subsequent Xn = ( ) 

points, as magenta box symbols, of the iteration sequence x,4; = f (Xn), : 

where f(x) = Ax. A solid magenta box distinguishes x,y, the final point for n = 0,1,2,.... 

calculated. Mathcad also calculates and plots the eigenlines of the 

matrix A (if there are any); if there are two eigenlines and one is 

dominant, then the dominant one is shown in red, and the other in blue. 

You can rescale the graph by altering the value of s. 


You can set A to be a matrix representing one of the basic linear In this activity, you will only 
transformations (or a product of such matrices). Alternatively, you can set A to be a scaling. 
set A to be the 2 x 2 matrix M, whose entries you can edit. 


The page is set up with A = S(1.5,1), (ao, yo) = (2,1), N=lands=25. Figure 5.1 illustrates the 

(a) Set N to 2, 3 and 4 in turn, so the third, fourth and fifth points of the see = aes by this 
sequence appear on the graph. (Making the points appear one at a Sees & amr Su ta 
time shows the order in which they appear in the sequence.) 

(b) Set N to 10, 20 and 30 in turn, and describe the effect on the value of | When you do this, only a few 
the ratio yy /xy, which is displayed to the left of the graph. Explain points of the sequence appear 
how you could have predicted this. on the graph, because the 


ee — : b t points li tsid 
Set the initial point (xo, yo) to (—2,1). Observe that the iteration ree ete eae eS es 


the graph box. You can 
sequence now tends in the negative direction of the z-axis. Se ae ee 


“——™ 
O 
Stn” 


Experiment with a few other initial points of your own. In each case, points, if you wish. 
predict the behaviour of the sequence and check your prediction. 


Reset the initial point to (2,1), and ensure that N = 30 and s = 25. 


Poti 
0 
wae 


Set A to each of the non-uniform scaling matrices below in turn, and 
in each case check that the behaviour of the sequence is as predicted 
by the iteration properties of generalised scalings. For the last of these matrices, 


] . 
S(1.5, 1.2), S(1.5, 0.8), S(1.5, = 12), SUS, —0.8), S(0.8, —0.8). neg the graph by setting 


Solutions to parts (b) and (c) are given on page 35. 


2 


See Chapter B3, Section 3. 
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Mathcad notes 


The solid box that marks the point (xy, yy) on the graph is obtained by 
plotting three traces, using a box, a ‘x’ and a ‘+’. (The traces are plotted 
using weight ‘p’, which draws the smallest symbols possible.) When 
superimposed on the screen, these three symbols give the appearance of a 
solid box. The individual symbols may become visible if printed. 


In the next activity you will consider iteration sequences produced by 
generalised scalings for which the eigenlines are different from the axes. 
You have already seen some examples of such sequences in the main text. 
For example, you considered the linear transformation with matrix 


i 2 
= = = 
3 


which has eigenvalues 4 and —1, with corresponding eigenlines y = 5x and 
y = —«, respectively. By the iteration properties of generalised scalings, 
the long-term behaviour of iteration sequences whose initial point is not on 
the eigenlines can be described in this case as follows. 


The points of such an iteration sequence alternate between opposite sides 
of the eigenline y = 2 and lie on the same side of the eigenline y = —z as 
the initial point. The sequence moves away from (0,0) and tends in the 


direction of the dominant eigenline y = $a. 


Figure 5.2 shows the first few points of the iteration sequence with initial 
point (2,1). 


Figure 5.2 Iteration sequence generated by a generalised scaling 


The next activity involves examples of a similar kind. 


SECTION 5 Iterating linear transformations with the computer 


Activity 5.2. Linear transformations with two distinct 


eigenvalues 
Ensure that the number of iterations, initial point and graph scale on You should still be working 
page 3 of the worksheet are set as follows: with Mathcad file 221B3-01. 


N = 30, eo Fe! s = 100. 
Yo 1 


In each of parts (a) to (d) below, a matrix A is given, together with its 
eigenvalues and eigenlines. In each case, use this information together with 
the iteration properties of generalised scalings to predict the behaviour of 
the iteration sequence x, = A”X with initial point (2,1), giving a 
description along the lines of the example discussed before this activity. 
Then confirm your answer by setting A appropriately in the Mathcad page 
and observing the graph. 


In part (a), define A := M; the matrix M is already entered for you. For 


the matrices in parts (b) to (d), edit the entries of M. For some of these matrices, 
Se you may find it helpful to set 
(a) ( — ) N to 1, 2, 3 and so on, to see 
’ the order in which the points 
has eigenvalues 2 and —3, with eigenlines y = 4x and y = —az, aitaiies 
respectively. sn 
i 
0 (23) 
has eigenvalues 4 and 2, with eigenlines y = x and y = si, respectively. 
(c) ( oe “e For parts (c) and (d), set 
: : = 3. 
has eigenvalues 0.7 and 0.6, with eigenlines y = —2z and y = —3z, : 
respectively. 
—0.6 0.5 
(d) ( 0.9 0.6 ) 
has eigenvalues 0.9 and —0.9, with eigenlines y = 3x and y = —2r, 
respectively. 


Solutions are given on page 35. 


In the next activity you will see some examples of sequences obtained by 
iterating linear transformations whose matrices have no eigenvalues. 
(Recall that when we say in this course that a matrix has no eigenvalues, 
we mean that it has no real eigenvalues. It will have complex eigenvalues, 
but we are not concerned with those.) Such sequences show a wide variety 
of behaviour, and you are not expected to explore them in detail. 
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You should still be working 
with Mathcad file 221B3-01. 


For some of these matrices, 
you may find it helpful to set 
N to 1, 2, 3 and so on, to see 
the order in which the points 
appear. 


When setting A to a product 
of two matrices, remember to 
type * (or use the ‘x’ button 
on the ‘Calculator’ toolbar) 


to indicate the multiplication. 
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Activity 5.3 Linear transformations with no eigenvalues 


Ensure that the number of iterations, initial point and graph scale on 
page 3 of the worksheet are set as follows: 


N = 200, (= 3=(5): eae 
Yo 1 


The matrices listed in parts (a) to (d) below have no eigenvalues; nor 
therefore do they have any eigenvectors. For each matrix, use your 
geometric knowledge of basic linear transformations to try to predict the 
general ‘shape’ of the iteration sequence produced by the matrix. (This is 
not easy for the matrices in parts (c) and (d)!) Then confirm your answer 
by setting A appropriately in the Mathcad page. 


(a) Rotations: 
RG), RG), RQ). 


6 3 
(b) Rotations composed with uniform scalings: 
U(0.9) R(4), U(-0.9)R(8), U(1.01) R(¥)- 
(c) A rotation composed with a non-uniform scaling: 
$(0.8, 1.2) R(4). 
(d) Rotations composed with shears: 
X(0.2)R(S), X(-0.6) R(B). 
Solutions are given on page 36. 


Comment 


No matrix of the types in parts (a) and (b), representing rotations, or 
composites of rotations and uniform scalings, has any eigenvalues (except 
where the angle of rotation is an integer multiple of 7). However, some 
matrices of the types in parts (c) and (d), representing composites of 
rotations and non-uniform scalings, or composites of rotations and shears, 
can have two eigenvalues. For example, $(0.8, 1.4) R(4) and X(0.3) R(,5) 
both have two eigenvalues. 


In the final, optional, activity in this subsection, you can see a few 
examples of sequences obtained by iterating linear transformations whose 
matrices have just one eigenvalue. Again, such sequences show a wide 
variety of behaviour, but it is worth looking at one particular aspect of this 
behaviour. You have seen that if A is a matrix with two eigenvalues of 
different magnitudes, and the initial point x9 does not lie on an eigenline 
of A, then 


=~ + mi 4b M+ 60, 


Ln 
where (pn, Yn) is the (n + 1)th point of the sequence x, = A”X and m is 
the gradient of the dominant eigenline. In the next activity you are invited 
to explore whether a similar property holds for 2 x 2 matrices with just 
one eigenvalue. 


SECTION 5 Iterating linear transformations with the computer 


Activity 5.4 Linear transformations with one eigenvalue 
(Optional) 


Ensure that the initial point and graph scale on page 3 of the worksheet 
are set as follows: 


(2 )=(7). s=50 


Each of the six matrices in parts (a) and (b) below has just one eigenvalue 
and just one eigenline. In each case set N = 50, set A to the given matrix, 
and note the gradient of the eigenline. Then increase N to 500, and check 
whether the value of the ratio yy /xy seems to be approaching the gradient 
of the eigenline. 


For some of these matrices, you will need to rescale the graph to obtain a 
clear picture of the iteration sequence. 


(a) (i) X(4) (ii): U(1.5) X(4) (iti): U(0.5) X(4) 
a font 1 i at oud lb 2 
(b) @) = 3 my ( 1 oe — ( —0.2 <% 


Solutions are given on page 36. 


Comment 


Most 2 x 2 matrices with just one eigenvalue have just one eigenline. The 
only exceptions are the matrices of uniform scalings and the zero matrix; 
for these matrices, every line through the origin is an eigenline. 


It can be shown that every linear transformation of the plane whose 
matrix has only one eigenvalue is a composite of a shear and a uniform 
scaling. If the matrix has only one eigenline, then the shear is parallel to 
this eigenline. 


Mathcad notes 


In part (a)(iii), you may have noticed that for both N = 50 and N = 500 
Mathcad displays the last point calculated as (ay, yn) = (0,0) but gives a 
non-zero value for the ratio yy /xy, despite the fact that Mathcad 
(somewhat dubiously) evaluates 0/0 as 0. The reason for this is that, while 
Mathcad always retains values to 15 significant figures for calculation 
purposes, in this worksheet numbers of magnitude less than 10~'° are 
displayed as zero on the screen. So Mathcad displays the values of 

59 = 1.794 x 107*8 and ys9 = 8.882 x 107° as zero, while displaying the 
calculated value of ys0/a59 as 4.95 x 107°. 


If you wish to display such small values, then click in an empty space in 
the worksheet to obtain the red cross cursor, and select Result... from the 
Format menu. In the ‘Result Format’ option box, choose the ‘Tolerance’ 
tab, and increase the value of ‘Zero threshold’. For example, increasing 
this value from 10 to 100 will ensure that only numbers of magnitude less 
than 10~'°° are displayed as zero. 


Now close Mathcad file 221B3-01. 


You should still be working 
with Mathcad file 221B3-01. 


The gradient of the eigenline 
is the y-component of the 
eigenvector with z-component 
equal to 1. 


In part (b), define A := M 
and edit the entries of M. 


One or both of the 
transformations forming the 
composite may be the 
identity transformation. 
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If you are not sure about the 
order in which the points 
appear in a sequence, then set 
N to 1, 2, 3 and so on, in 
turn, to find out. 
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5.2 Iterating affine transformations 


In Subsection 5.1 we looked at sequences of points in the plane obtained by 
iterating linear transformations. In the next activity you will see a few 
examples of iteration sequences obtained by iterating affine 
transformations. You saw in Chapter B2 that an affine transformation of 
the plane is a function of the form 

f:R? — R’ 

xt Ax-+a, 

where A is a 2 X 2 matrix and a is a vector with two components. | 


We generate an iteration sequence using an affine transformation f in a 
similar way to the other iteration sequences that you have seen: we choose 
an initial point x») and repeatedly use the recurrence relation 


4.11 = Pie (e842; 


Activity 5.5 Affine transformations 


Open Mathcad file 221B3-02 Iterating affine transformations. Page 2 
of the worksheet defines the matrices representing basic linear 
transformations, using the Mathcad notation. 


Move to page 3. This is similar to page 3 of the worksheet for Mathcad 
file 221B3-01, but it is set up to allow you to define a vector a as well as a 
matrix A, and the function iterated is f(x) = Ax +a. The initial point 
(xo, Yo) is set to (2,1), and the number N of iterations is set to 10. 


(a) The matrix A is set to the rotation matrix R(7), and the vector a is 


0 ee ; 
set to ( 0 ), so f is initially a linear transformation. 


Set N = 20 and describe the effect on the iteration sequence. 


Then set a to each of the following vectors in turn and describe the 
effect on the iteration sequence: 


0 2 a 
| Ee & —1 }) 
(b) Set A = U(0.9) R(2), which represents a rotation composed with a 


uniform scaling, and set a = ( : ). 


Set N = 40 and describe the effect on the iteration sequence. 


Then set a to each of the vectors in part (a) and describe the effect on 
the iteration sequence. 


Solutions are given on page 36. 


Now close Mathcad file 221B3-02. 


SECTION 5 Iterating linear transformations with the computer 


It can be shown that if the 2 x 2 matrix A represents an anticlockwise 

rotation through an angle @ about the origin, then each affine 

transformation f(x) = Ax +a is an anticlockwise rotation through the 

angle 6 about some point whose coordinates depend on @ and a. Thus if A’ You saw in Chapter B2, 
represents a rotation, then each affine transformation f(x) = Ax+a Section 4, how to find a so 


produces iteration sequences whose points lie on circles. that f is a rotation about a 
ges given point. 
Similarly, if A represents a composite of a rotation and a uniform scaling, 


then each affine transformation f(x) = Ax + a produces iteration 
sequences that spiral about a point whose coordinates depend on @ and a. 


5.3 Iterated function systems (Optional) 


In Subsection 5.2, you looked at some sequences of points in the plane Do not be tempted to spend 
obtained by iterating affine transformations. In this subsection, you will too long on this subsection! 
see some sequences obtained in a similar way, but using more than one 

affine transformation. A collection of transformations used to generate an 

iteration sequence in this way is known as an iterated function system. 

Such systems can give surprising results! See Activity 5.7, for example. 


For example, consider the two affine transformations 


f(x) =Ax+a and g(x) = Ax-a, 


where A = U(0.6) R(0.5) and a= ( oe ) 


0 
The matrix A represents a composite of a rotation and a uniform scaling. The angle of rotation is 
Both f and g produce spiral sequences if iterated on their own, the spirals 0.5 radians (about 29°) and 
centred on points other than the origin, since a is not the position vector the scaling factor is 0.6. 


of the origin. The sequences spiral inwards towards their centres, because 
the scaling factor has magnitude less than 1. 


We can produce an iteration sequence using both f and g by starting with 
an initial point, and choosing in some way between f and g at each 
iteration. For example, we could simply choose either f or g at random. 
Figure 5.3 shows a plot of the first 50000 points of a sequence obtained by 
iterating f and g in this way; the initial point is the origin. 


“4 “0.5 0 0.5 l 


Figure 5.3 Sequence generated by an iterated function system 


oF 


Complicated self-similar sets 
are often called fractals. 


In part (b), Mathcad starts a 
new plot after each change. 
For ways to avoid waiting, see 
the Comment on page 13. 


32 


CHAPTER B3 


As you can see, the sequence produces a plot of a set that is of great 
complexity yet displays the property of self-similarity; that is, small parts 
of the set are ‘similar’ to large parts. Moreover, the set has an interesting 
‘natural’ appearance, perhaps resembling some kind of sea life! 


In the next activity you are invited to use a Mathcad page to see a variety 
of plots obtained in a similar way. The page allows you to define two affine 
transformations f and g, and an initial point x9. It generates and plots a 
sequence x, using the following method. At each iteration, Mathcad 
generates a random number p,, between 0 and 1, and the next point of the 
sequence is calculated using the recurrence relation 


= fie Bee 
whl 7 | ga i pee 


where P is a number defined in the worksheet. The variable P is initially 
set to 0.5, so each of f and g is chosen approximately equally often. 


Activity 5.6 Two-function systems (Optional) 


Open Mathcad file 221B3-03 Iterated function systems. Page 2 of the 
worksheet defines the matrices representing basic linear transformations, 
using the Mathcad notation. 


Move to page 3. Here two affine transformations 
f(x) =Ax-+a and g(x) =Bx+b, 

are defined, where 
A =U(r1) R(61) and B=U(r2) R(62). 


The variables r1, 61, r2 and 02, and the vectors a and b, are initially set 
to values that give the transformations f and g used to generate the 
sequence in Figure 5.3. 


The initial point xo is set to the origin. The iteration sequence x,, of 
points is generated by using a sequence p,, of random numbers and a 
variable P to decide which of the transformations f and g is used at each 
iteration, in the way described before this activity. The number of 
iterations carried out is NV. 


Each calculated point of the sequence is displayed on a graph as a small 
dot. The graph can be rescaled by changing the values of s1 and s2. 


The number of iterations, N, is set to 5000. If your computer is fast 
enough, then you may wish to increase this to 20000 or 50000, for 
example, to obtain ‘better’ plots. 


(a) Make the following three changes in turn (do not make other changes 
as’you do this), and observe each new plot obtained. 


(i) rl1=0.85 (ii) rl1=0.7 (iii) 01 =0.3 
(b) Make the following changes, and observe the plot obtained after they 
are all made: 


rl=095, #1=05, 2=08, P=09. 
(You may also wish to set s1 = —1.) 


Comment 


By varying the parameters r1, 01, r2 and 62, you can alter the shape of 
the sequence generated by the iterated function system and so produce a 
wide variety of ‘natural-looking objects’. 
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Mathcad notes 


© The Mathcad function if chooses one of two values depending on a 
condition. The expression if (condition, a,b) gives the value a if the 
condition is true, and 6 if the condition is false. Several if statements 
can be combined together to make multiple choices. 


© The Mathcad function rnd generates uniformly-distributed random 
numbers. The expression rnd(z) returns a random number between 0 
and x. The values produced by the rnd function come from a 
sequence of random numbers generated by Mathcad when it starts up. 
The same sequence is generated every time, unless you change the 
random number generator. To do this, select Options... from the 
Math menu, choose the ‘Built-In Variables’ tab and then set the ‘Seed 
value for random numbers’. Each choice of positive integer gives a 
different sequence; the default value is 1. 


© ‘The graphs in the worksheet are plotted using the trace type ‘points’, 
with symbol ‘none’. Mathcad plots each point as a small dot. 


It is a remarkable fact that the shape of the set obtained by this method is 
associated only with f and g, and not with the way that we choose 
between them at each iteration. Briefly, if the functions f and g are both 
contracting, that is, they reduce the distance between pairs of points, then 
there is a bounded set F of points which has the property that if we 
choose any point in F,, and apply either f or g to this point, then we 
obtain another point in F’. If we choose a point outside F’, and repeatedly 
apply either f or g, then the resulting sequence approaches the set F', and 
may or may not eventually lie within F’. 


One way to obtain a picture of the set F' is to iterate f and g in the 
random manner described. This usually produces a ‘chaotic’ sequence 
whose points, at least after the first few, are either in the set or very close 
to the set. For some pairs of transformations f and g, such as the pair in 
Activity 5.6(b), we need to ensure that one of the transformations is 
applied more often than the other if we want the points of the sequence to 
be fairly evenly distributed across the set. 


The self-similar nature of sets like those in Activity 5.6 led mathematicians 
to the remarkable discovery that it is possible to use iterated function 
systems to generate a given set by analysing the ways in which it is 
self-similar. This usually requires more than two affine transformations. A 
classic example is the fern, which you can see in the final activity. 


Activity 5.7 The fern: a four-function system (Optional) 


Move to page 4 of the worksheet, and observe the plot obtained from the 
given iterated function system of four affine transformations. 


If your computer is fast enough, then you may wish to increase the number 
of iterations N to 20000 or 50000, for example, to obtain a better plot. 


If you wish to experiment with the plot, then you could change the 
rotation matrix in the definition of the matrix A; for example, try R(0) 
and then R(0.05). Other small changes to the parameters in the matrices 
will give ferns of slightly different shapes. 


Now close Mathcad file 221B38-038. 


See A Guide to Mathcad. 


You should still be working 
with Mathcad file 221B3-03. 
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Solution 4.2 


(a) When c takes any of the values 0.5, 0.6 or 0.7, 
the sequence tends to infinity. 


(b) When c takes either of the values —0.1 or 0.1, 
the sequence tends to a fixed point. 


(c) When c takes any of the values —0.5, —0.6 or 
—0.7, the sequence tends to a fixed point. 


(d) When c takes the value —1.3, the sequence tends 
to a 4-cycle. 
When c takes either of the values —1.1 or —1.2, 


the sequence tends to a 2-cycle. 


(e) When c takes any of the values —1.7, —1.8 
or —1.9, the sequence seems to be chaotic. 


The results of parts (a) to (e) are summarised in the 
following figure. 


5 4-2-3 ae ee 
—1.9 —1.7 —1.3 —1.1 —0.7 —0.5 
—1.8 —1.2 —O0 
FF 00 26 00 


Solution 4.3 


(a) The sequence z,, tends to a fixed point for all 
values of c in the interval [—0.75, 0.25]. 


(For values of c inside this interval but close to 
its left end, this convergence is difficult to see 
from the graphical iteration diagram or the 
table, but it is confirmed by the fact that the 
gradient of the relevant fixed point lies between 
—1 and 1, so the fixed point is attracting. Note 
that if c is inside the interval but very close to 
either end, then the gradient at the attracting 
fixed point is displayed as +1, due to rounding. 


When c = 0.25 and c = —0.75 the relevant fixed 
points are indifferent, and in both cases you 
have to set Mathcad to carry out a large number 
of iterations to see that the sequence appears to 
tend to the fixed point. 


For values of c a little greater than 0.25, the 
sequence tends to infinity. 


For values of c a little less than —0.75, both 
fixed points are repelling, and the sequence 
tends to a 2-cycle.) 
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(b) The sequence z,, tends to a 2-cycle for all values 
of c in the interval [—1.25, —0.75). 


(When c = —1.25 the relevant 2-cycle is 
indifferent, and you have to set Mathcad to 
carry out a large number of iterations to see that 
the sequence appears to tend to the 2-cycle. For 
values of c a little less than —1.25, the 2-cycle is 
repelling, and the sequence tends to a 4-cycle.) 


Solution 4.4 


(a) The behaviour of the sequence z,, for the given 
values of c is as follows: 


—1.39: 2x, tends to an 8-cycle; 
—1.475: 2x, tends to a 6-cycle; 
—1.575: 2x, tends to a 7-cycle; 
—1.76: 2x, tends to a 3-cycle. 


(b) Setting c = —1.4 produces a graphical iteration 
diagram with a fairly large number of 
construction lines. Increasing the number of 
iterations does not produce any visible new 
construction lines, so it appears that the 
sequence tends to a cycle. The number of 
members of the cycle cannot be counted in the 
diagram scaled from —2 to 2, since some of the 
construction lines are too close together. (In 
fact, this is a 32-cycle.) 


Solution 4.6 


(b) Values of c between about —1.628 and —1.625 
give 5-cycles. 


Solution 4.7 


(a) The three fixed points are 0, 0.766 and —0.766 
(to three decimal places). They are all repelling. 


(b) The two 2-cycles are —1.034, —0.467 and 
1.034, 0.467 (to three decimal places). They are 
both attracting. 


(c) When 2p = —1, the sequence tends to the first 
2-cycle given in the answer to part (b). 


When zo = 1, the sequence tends to the second 
2-cycle given in the answer to part (b). 


Chapter B2 


Solution 5.2 


(c) (i) The matrix represents the composite of a 
scaling with factors 2 and 1, followed by an 
x-shear with factor 2; that is, 


( an ~ X(2) $(2,1) 


(You can check the answer by matrix 
multiplication. ) 


(ii) The matrix represents the composite of a 
scaling with factors 2 and 1, followed by an 
x-shear with factor 2, followed by a rotation 
through the angle 7; that is, 


( % a = R(m) X (2) $(2, 1) 


=(% 4)(o 1) (0 1): 


(iii) The matrix represents the composite of an 
x-shear with factor —1, followed by a rotation 
through the angle 7/2; that is, 


(ft) =R@) xen 


=(1 o)(o 1) 


(iv) The matrix represents the composite of a 
scaling with factors —3 and 3, followed by an 


x-shear with factor —3, followed by a rotation 


through the angle 22 (or, equivalently, — 5); that 


2 
is, 


Solution 5.3 
Target figure 1 is obtained by setting 


A=R(") and a=(}). 


Target figure 2 is obtained by setting 


A= QO) and a= ( : ). 
Target figure 3 is obtained by setting 


A. =3(2,3) and a= ( = . 


Target figure 4 is obtained by setting 


A=Y¥@)- and a={ 7 : 


SOLUTIONS TO ACTIVITIES 


Chapter B3 


Solution 5.1 


(b) 


The ratio yy/xy appears to tend to 0 as N 
increases. ‘This could have been predicted 
because 0 is the gradient of the dominant 
eigenline y = 0. 


If the initial point lies to the right of the y-axis, 
then the sequence tends in the positive direction 
of the x-axis, whereas if the initial point lies to 
the left of the y-axis, then the sequence tends in 
the negative direction of the z-axis. 


(The dominant eigenvalue 1.5 is positive, so each 
point of the sequence lies on the same side of the 
non-dominant eigenline as the initial point.) 


If the initial point lies on the y-axis, then 
subsequent points of the sequence coincide with 
the initial point. (In this case, Mathcad 
highlights the ratio yy /xy in red and displays 
no value for it, since its calculation would 
require division by zero.) 


Solution 5.2 


(a) 


The points of the sequence alternate between 
opposite sides of the eigenline y = 4x and stay 
on the same side of the eigenline y = —z as the 
initial point. The sequence moves away from 
(0,0) and tends in the direction of the dominant 
eigenline y = —2. 


The points of the sequence stay on the same side 
of the eigenline y = x and on the same side of 
the eigenline y = 2 as the initial point. The 
sequence moves away from (0,0) and tends in 
the direction of the dominant eigenline y = z. 


The points of the sequence stay on the same side 
of the eigenline y = —2z and on the same side of 
the eigenline y = —3z as the initial point. The 
sequence moves towards (0,0) and tends in the 
direction of the dominant eigenline y = —2z. 


The points of the sequence stay on the same side 
of the eigenline y = —3r as the initial point, and 
alternate between opposite sides of the other 
eigenline y = 3x. The sequence moves towards 


(0,0). Neither eigenline is dominant. 


(In fact, the points lie alternately on each of a 
pair of lines through the origin.) 


a5 


COMPUTER BOOK B 


Solution 5.3 


(a) The matrix R(%) gives a sequence each of whose 
points is one of twelve points equally spaced 
around a circle centred at the origin. 


The matrix R(%) gives a sequence each of whose 
points is one of six points equally spaced around 
a circle centred at the origin. 


The matrix R(1) gives a sequence whose points 
lie around a circle centred at the origin; no 
points are repeated. 


(b) The matrix U(0.9) R(;5) gives a sequence whose 
points lie on a spiral centred at the origin; each 
point is closer to the origin than its predecessor. 


The matrix U(—0.9) R(#5) gives a sequence 
whose points lie alternately on each of two 
spirals centred at the origin; each point is closer 
to the origin than its predecessor. 


The matrix U(1.01) R(75) gives a sequence 
whose points lie on a spiral centred at the origin; 
each point is further from the origin than its 
predecessor. 


(c) The matrix 5(0.8, 1.2) R(75) gives a sequence 
whose points lie on an ‘elliptical spiral’, centred 
at the origin and spiralling towards the origin. 


(d) The matrix X (0.2) R() gives a sequence whose 
points lie on an ellipse centred at the origin. 


The matrix X(—0.6) R(75) also gives a sequence 
whose points lie on an ellipse centred at the 
origin. 


Solution 5.4 


The value of the ratio yy /xx seems to approach the 
gradient of the eigenline in all six cases. 


(It is indeed true that if A is a 2 x 2 matrix with just 
one eigenvalue and just one eigenline, and Xo is a 
point in the plane, then 


oo in ae OS 


nr 
where (2p, Yn) is the (n+ 1)th point of the iteration 
sequence x, = A”x and m is the gradient of the 
eigenline. ) 
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Solution 5.5 


(a) 


Setting N = 20 illustrates that after the first 
twelve points the sequence begins to repeat. 
Each point in the sequence is one of twelve 
points equally spaced on a circle centred at the 
origin. (This sequence was the first of those 
considered in Activity 5.3(a).) 


Changing a to a vector other than the position 
vector of the origin produces a sequence whose 
points lie on a circle centred at a point other 
than the origin. 


(Changing a usually also changes the radius of 
the circle, since this is the distance between the 
centre of the circle and the initial point.) 


Setting N = 40 plots more points of the 
sequence, spiralling towards the origin. 


Changing a to a vector other than the position 
vector of the origin produces a sequence that 
spirals inwards about a point other than the 
origin. 
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